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Abstract 

Caratheodory showed that n complex numbers c±, . . . , c n can uniquely 
be written in the form c p = 2j=i Pj e j p with p = 1, . . . , n, where the 
€jS are different unimodular complex numbers, the pjS are strictly pos- 
itive numbers and integer m never exceeds n. We give the conditions 
to be obeyed for the former property to hold true if the pjS are sim- 
ply required to be real and different from zero. It turns out that the 
number of the possible choices of the signs of the pjS are at most equal 
to the number of the different eigenvalues of the Hermitian Toeplitz 
matrix whose i, j-th. entry is Cj-i, where c_ p is equal to the complex 
conjugate of c p and cq = 0. This generalization is relevant for neutron 
scattering. Its proof is made possible by a lemma - which is an inter- 
esting side result - that establishes a necessary and sufficient condition 
for the unimodularity of the roots of a polynomial based only on the 
polynomial coefficients. 
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1 Introduction 



Caratheodory's theorem 1 states that n complex number Ci, . . . , c n as well as 
their complex conjugates, respectively denoted by c_i, . . . , c_ n , can always 

and uniquely be written as 

in 

c P = J2pj e /i P = 0,±l,...,±n, (1) 

i=i 

with pj G R, pj > 0; 6j G C, \ej\ = 1, €j ^ €k (j ^ k = 1, . . . , m) and c 
uniquely determined by the CjS with j 7^ 0. Further, the pjs, the e^s and m 
are unique and m obeys to 1 < m < n. 

The practical relevance of this theorem for the inverse scattering problem 2 as 
well as for information theory 3 ' 4 appears evident from the following remark. 
Writing the e^s as e l2nx ^> with < Xj < 1, the c p 's take the form J2]Li Pje t2lTXjP 
so as they can be interpreted as the scattering amplitudes generated by m 
point scatterers (with "charges" p±, . . . , p m respectively located at xi, . . . , x m ) 
and relevant to the "scattering vector" values p — 0, ±1, . . . , ±n. One con- 
cludes that the theorem of Caratheodory allows us to determine the positions 
and the charges from the knowledge of the scattering amplitudes Ci, . . . , c„. 
Further, it ensures that the solution of this inverse problem is unique. How- 
ever, in the case of neutron scattering 5 , the charges of the scattering centers 
have no longer the same sign. Nonetheless, it has recently been shown 6 that 
the so-called algebraic approach for solving the structure in the case of X-ray 
scattering from an ideal crystal 2 can also be applied to the case of neutron 
scattering. This result suggests that the aforesaid Caratheodory theorem 
can be generalized so as to avoid the requirement that the sign of all the 

2 



scattering charges be positive. In this note we show how this generalization 
is carried out. 

Before proving this statement, we find it convenient to briefly sketch the proof 
of Caratheodory's theorem reported by Grenander and Szego 7 . The proof is 
based on an enlargement of the set of the n given complex numbers c p to a 
set containing (2n + l) complex numbers still denoted as c p with — n < p < n, 
the c p s with negative index being defined as the complex conjugates of the 
given c p s, i.e. c_ p = Cp~ with p = 1, . . . , n. (Hereafter an overbar will always 
denote the complex conjugate). The remaining value Cq, real by assumption, 
is determined as follows. Consider the (n + 1) x (n + 1) matrix (C) with its 
(r, s)th element defined as 



This matrix is a Hermitian Toeplitz matrix 8 and its diagonal elements are 
equal to cq. This value is chosen in such a way that the matrix (C) turns out 
to be singular (i.e. det(C) = 0) and the associated bilinear Hermitian form 

n+l 

u\C)u= UrC rjS U s (3) 

r,s=l 

[where u is an (n + 1) dimensional complex vector] non- negative (or semi- 
positive) definite. To show that c can uniquely be determined, one proceeds 
as follows 3 . Consider the (n+l) x (n+l) matrix (C) that has its non-diagonal 
elements equal to the correspondent elements of (C) and the diagonal ele- 
ments equal to zero, i.e., with r, s — 1, . . . , (n + 1), 




(2) 
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This matrix is Hermitian. Then its eigenvalues (xj, j = 1, . . .n + 1), are 
real and can be labeled so as to have Xi ^ • • • ^ X n +i ■ Further, they are 
such that YZ[=i X 3 — because the trace of (C) is zero. Hence, Xi < an d 
X n+ i > 0. One immediately realizes that matrix (C) is obtained by setting 
c = — Xi > so that (C) = (C) — Xi(-0j (-0 being the unit matrix. In 
fact, the matrix (C — X\I) is a Toeplitz Hermitian matrix with its diagonal 
elements equal to (— xj. The secular equation of this matrix is 

det((C) - z{I)) = det(((7) - - z(I)) = det((C) - ( Xl + *)(!)) = 0. 

This equation is the same equation that determines the eigenvalues of (C) if, 
instead of variable z, we use the shifted variable z + x 1 - Thus, the eigenvalues 
of (C) are: = (Xi - Xi) < (x 2 ~ Xi) < • • • < (x„+i ~ xj an d th e matrix 
(C) is semi-definite positive. Let /i x (> 1) denote the multiplicity of the 
eigenvalue x± °f then the rank of (C) is (n + 1 — fix) and the m present 
in Eq. (0) has the same value, i.e. m = (n + 1 — /ij. Exploiting the non- 
negative definiteness of (C), Grenander and Szego showed that: I) the e,-s 
are distinct and are the roots of the following polynomial equation, referred 
to as resolvent equation in the following, 



/ 



P m (z) = D m - 1 det 



■"0 



^-m+1 

1 2 



V 



\ 



„m— 1 



(5) 



where D m > is the determinant of the left principal minor contained in the 
first m rows; II) the pjS are strictly positive and given by 

i m— 1 



P.! 



P 'fe , 



(6) 



where the prime denotes the derivative and the f3j }P s are the coefficients of 
the following polynomial 

m—l 

P m (z)/(z - ej) = £ f3 hpZ v- (7) 

and III) that Eq. ((TJ holds true with p = —n, (—n + 1), . . . , n. 
This theorem, via implies that any non- negative definite Hermitian 
Toeplitz matrix (C), defined by (J2J), uniquely factorizes as (V)(Q)(V)^ where 
(V) is an (n + 1) x m Vandermonde matrix with Vjj = , (V)^ its Her- 
mitian conjugate and (Q) an m x m positive-definite diagonal matrix with 
Qi j = piSij. A more general factorization was obtained by Ellis and Lay 9 . 
In their theorem 3.4, they showed that a Toeplitz matrix of order (n + 1) and 
rank m <n can be factorized as (V)(A)(V)^ where (V) is a generalized con- 
fluent Vandermonde matrix generated by the distinct roots of the resolvent 
equation and (A) is a block- diagonal matrix with reversed upper triangular 
blocks. In proving this generalized factorization, the Toeplitz matrix was 
neither required to be Hermitian nor positive/negative semidefinite and this 
last property - as it will appear clear in § 2 - amounts to relax the condition 
that the pjS present in (JJ) be strictly positive. However, in proving this the- 
orem, it was explicitly assumed that all the roots of the resolvent equation 
were unimodular without stating the conditions that have to be obeyed by 
the Toeplitz matrix elements for the unimodularity property to occur. On 
the one hand, this point makes the generalization not fully proven. On the 
other hand, the physical motivation that underlies our attempt of general- 
izing Caratheodory's theorem requires that Toeplitz matrices are Hermitian 
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as well as a strictly diagonal structure of (A). 

The generalization of Caratheodory's theorem presented below meets the last 
two requirements. In particular, in appendix A we report the lemma that 
specifies the necessary and sufficient conditions that must be obeyed by the 
coefficients of a polynomial equation for all its roots to lie on the unit circle. 
We also stress that all the so far reported 10,11 theorems, ensuring the unimod- 
ularity of the roots, are not based on the only knowledge of the polynomial 
coefficients, as our lemma does. Finally, appendices B, C and D report some 
properties of the coefficients of the resolvent equation, of Hermitian Toeplitz 
matrices and some numerical illustrations, respectively. 

2 Generalization of Caratheodory's theorem 

The aforesaid choice of Co is the only one that yields a non negative definite 
matrix (C). It is evident that if we choose c = — x n+1 the resulting (C) is 
a non-positive definite matrix and — (C) is a semi-positive definite one. We 
can apply to the latter the same analysis made by Grenander and Szego with 
the conclusion that 

m 

-° P = p'j e 'j p i p = o, ±1, . . . , ±n 

i=i 

where m — (n+1 — /i n+1 ) , (i n+1 denoting the multiplicity of x n +i ? the p' s are 
positive and the e'js are distinct unimodular numbers. 

Of course this generalization is quite trivial. The remaining choices are how- 
ever more interesting. In fact, if we choose Cq = — Xj with Xj 7^ Xi an d 
Xj 7^ Xn+i j the resulting matrix (C) will be indefinite, since some of its eigen- 
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values are negative, others positive and one, at least, equal to zero. The 
question is now the following: does a relation similar to (JJ hold also true in 
this case? We shall show that the answer is affirmative iff: i) m, the rank of 
(C) smaller than the order (n + 1) of (C), is equal to the principal rank of 
(C) and ii) the conditions reported in the lemma discussed in appendix A are 
fulfilled. In i) principal rank means the highest of the rank values of all the 
strictly principal left minors of (C) where, by definition, an (m x m) strictly 
principal (left) minor of (C) is a principal minor formed by m subsequent 
rows and columns with indices p + 1, . . . , p + m and < p < (n + 1 — m). 
Let x 1 < . . . < \ v denote the distinct eigenvalues of (C) (the matrix with 
Co = 0), fi 1 ,...,fi h , their respective multiplicities (independent from the 
choice of Co) and and Py. the number of the xis that respectively are 
negative and positive, so that either v = (n% + p^ + 1) or v = {n^ + p^) 
depending on whether one of the Xi s is equal to zero or not. If we set 

Q,o = ~Xh Ni_ = ^2^ and = J2 (8) 

t=l i=l+l 

the resulting matrix (Ci) will have Ni_ negative eigenvalues, iVj )+ positive 
eigenvalues and fii eigenvalues equal to zero so that its rank is mi = {n + 1 — 
/iz). We shall prove that the following relations 

mi 

c p = PiJ e i/i p = ±1, • • ■ , ±n, (9) 

3=1 
mi 

Q,o = J2pij, P = 

i=i 

hold true if mi is equal to the principal rank of (C/) and the conditions 
stated in the lemma are obeyed. In (jH)) the value of mi has been just defined, 
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the 6ijS are the distinct unimodular roots of a polynomial equation of degree 
mi uniquely determined by the c p s and c^o and, finally, the number of the 
positive and negative pijs respectively is Ni )+ and Ni-. 
It is remarked that we can construct v different matrices (Q) from a given 
(C) by setting q i0 = — Xi, I — 1, . . ■ ,v. Among these matrices only those 
which obey i) and ii) allow us to write Q 5 o and the c p s, with p — ±1, . . . , ±n, 
in the form (jHJ). Hence, if the number of these matrices is denoted by /i, one 
concludes that we have /i different ways for writing the n complex numbers 
(\, . . . , c n in the form 0. In general [i is such that 2 < \i < v. Clearly this 
result generalizes Caratheodory's theorem and we pass now to prove it. 
The proof must be achieved by a procedure different from that followed 
by Grenander and Szego because matrices (Ci) are no longer non- negative 
definite. Since our attention will focus on a particular (Ci), for notational 
simplicity we shall omit index / in the following considerations. Assume that 
Eqs Q are fulfilled so that 

m 

CV, S = c s - r = Pj e j~ r > r, s = 1, . . . , (n + 1), (10) 

with 

pj E R, pj ^ 0, €j e C, \ej\ = 1, €j ^ €i if j ^ i, j, i = 1, . . . , m. 

(11) 

Then C r ^ s = C Sjr and the Hermiticity of (C) is ensured. Introduce now two 
further matrices: (V) and (Q) such that 

Vj, r = e r f\ j = l,...,m, r = 1, . . . , (n + 1), (12) 
Qj,i = PjSj,i, j, i = l,...,m. (13) 
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The first is an m x (n + 1) rectangular Vandermonde matrix and the second 
a diagonal mxm matrix. As the e^s are distinct, the rows of (V) are linearly 
independent. From the unimodularity of the 6jS and immediately follows 
that 

m m m 

J2 VK,iQi,jVj !S = Vi,rPAj V J,s = Y.Pi e l~ r = C r,s ( 14 ) 

i,j=l i,j=l i=l 

so that one has 

(C) = {V)KQ){V). (15) 

Consider now the p x p minor of (C) formed by the rows r\ < • ■ ■ < r p and 
the columns s± < • ■ ■ < s p . By Bezout's theorem 12,13 one finds that 

H c zzO- £ [n»M< v z::d ie i v z:i) (16) 

l<il<i2<---<i P <m j=l 

where, adopting Gantmacher's notation 12 , the lower and upper indices inside 
each determinant symbol denote the rows and the columns of the considered 
minors of (C), (V) and (V). This expression makes it clear that the deter- 
minant of any minor of (C) of order p > m is equal to zero because the order 
of (Q) is m 14 . For this reason the rank of (C) cannot exceed m. At the same 
time, if p = m, any (m x m) strictly principal minor det (C q+1 '" -' q+m (with 
0<q<n-\-l — p) of (C) will have determinant equal to 



[II Pi] II 



l<i<jr'<m 



(17) 



In fact, if p — m, the sum present in (|TT)j) involves a single term and, due to 
(1121) . in det(V' 1+1 • ' • q+rn ) we can factorize e\ in the first row, e\ in the second 

V 1,— ,m ' 

row and so on. The remaining matrix is a Vandermonde matrix so that 
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In the same way one shows that 

m 

d 4^r'r) = n^ n 

' ' '=1 l<*<i<m 

Finally, the unimodularity of the e^s yields Eq. (II 7)1 . In this way, it has 
been shown that, if the elements of (C) have the expression reported on the 
right hand side (rhs) of (|10|) and conditions are obeyed, matrix (C) has 
rank m and each of its principal minors formed by m subsequent rows is 
nonsingular. In this way the rank of (C) is equal to its principal rank and 
the necessity of condition i) is now clear 15 . Since these minors are equal, the 
first m rows (or columns) of (C) are linearly independent. One recalls that 
the polynomial equation having its roots equal to the m e,-s, can be written 
as 

m m 

P m (z) = 'R(z-e j ) = y E l a l J = (18) 

3=1 1=0 

with 

ai = (-r~ l E *i / = 0,...,(m-l) (19) 

l<ii<— <i m -i<m 

and a m = 1. From ([19)1 and the unimodularity of the e/S follows that a 
is unimodular and, therefore, is different from zero. From the condition 
Pm(^i) = follows that e™ = — Yjb=Q l o-i^i- After multiplying the latter by 
ej J , with j G Z, and setting q = m + j, one finds that 

m—1 

el = — E a i e i +9 ~ m > 5 G Z and i — 1, . . . , m. (20) 

1=0 

The substitution of these relations in (|rU|) yields 

m m—1 m—1 

c s - r = - E Pi w[ +s ~ T ~ m = - E aic s -r-{m-i)- (21) 

i=l 1=0 1=0 
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Taking (s — r) = 1,2, ... ,m, one obtains the following system of linear equa- 
tions 

c a + c l a l + c 2 a 2+ •• + c m-l 1m-l = — C m 

c_i a +c ai +ci a2+ •• + c m-2 Om-i — — c m _i ^22^ 

C-m+l Q-o +C- m +2 Ol +C- m +3 «2+ •• +Co a m -l = — Ci- 

This uniquely determines coefficients clq, . . . ,a m _\ of (fTHj) because the de- 
terminant of the coefficients in (|2*2*|) is different from zero and equal to D m , 
defined below Eq. (5). Recalling that a m = 1, resolvent equation (|18J) is fully 
determined and formally coincides with Grenander & Szego's Eq. (J3J). 
The solution of (|22j) is 

a { = (-) m ^ l detfC 1 '-- w+a '-' m+1 VD, l = 0,...,m (23) 

V I,--- ,m / 

D^det^C 1 - '• m ) ^0, (24) 

where in ()23j) we let I assume the value m since it yields reported 
below ()19j) . Further, setting / = in Eq. (|23|) and using (J58j) (with n = m + 1, 
p = 1 and q = 2), one finds 

a = (-) m det (C 2 - - m+1 ) /D = (-) m e^ . (25) 

V 1,- ,m / 

In this way one recovers that is unimodular as already reported below 
(HHJ). As discussed in appendix C, coefficients (J2H)) are partly related by the 
relations a7 = a m _j/ao, / = 0, ...,m. These imply that the roots of (|18j) 
obey to Ij~ = 1/ej with j = 1, . . . , m and i\, . . . , i m equal to a permutation 
of {l,...,m}. This result ensures the unimodularity only for those roots 
with ij = j. The unimodularity of all the roots is ensured by the further 
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conditions stated in the lemma of appendix A. Thus, the roots of (JTHJl are 
unimodular if the new (m + 1) x (m + 1) Toeplitz matrix (S), defined as 

S r)S = <j s - r , r, s = 1, . . . ,m + 1, (26) 

°p = Z e /'P = °' ±1 ' ±2 > • • • > ±m ' 

i=i 

is nonnegative definite and has rank m. Quantities a p are easily and uniquely 
determined by Eq.s (|39|) and (|40|) in terms of the a p s given by (|23|). Thus, 
after checking that the determinants of the left principal minors, contained 
in the first p rows of (S), are strictly positive for p = 1, . . . ,m and equal 
to zero for p = m + 1, the unimodularity of all the roots of (|18jl is ensured. 
After solving the resolvent equation, the pjS can be determined by Eq.s (jHJ 
and (J7J) since these also apply in the case of non-definite (C) . Alternatively, 
the pjS can be determined solving the system of m linear equations 

m 

£e/p, = c p , p = 0,...,(m-l) (27) 
i=i 

that follow from (JHJ). [For notational simplicity, we still omit index I present 
in the definition of Cq\. These equations can also be written as 

m 

Y. V P+ijPi = c p> P = 0,...,(m-1) (28) 

where {V T ) is the transpose of the mxm upper left principal minor of matrix 
{V) defined by flEJ). The formal solution of (J2EJ) is 

m— 1 

Pi = E (V T )- l Jp+l c p , j = l,...,m, (29) 

p=0 

since (V T ) is certainly non singular. 

Finally it must be proven that the numbers of the pjs that turn out to be 
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positive or negative are respectively equal to N + and N_ (again we omit 
index I). To this aim consider the Hermitian bilinear form 

n+l 

C sN = ^h-Cr^Us, u s e C. (30) 

r,s=l 

By Eq (|15jl this is immediately expressed in terms of the diagonal form 



C 2M = J2 v pi u ]pp v p[ u ] ( 31 ) 
P =i 



where 

n+l 



"pM = 51 ^>,s w s, p = 1, • • • , m. (32) 

8=1 

At the same time, since (C) is Hermitian, it can be diagonalized by a unitary 
transformation (U) and written as 

(C) = (l7)t(x)(E0 (33) 

where 

(X)r,s = (Xr ~ Xi)8r, s , r, s = 1, . . . , (u + 1) (34) 

with the XrS equal to the eigenvalues of (C). As discussed at the beginning 
of this section, fii of the (xr — Xi) s are equal to zero, A^ i+ are positive and 
Ni,- negative. Therefore we can compact (U>) by eliminating \xi columns 
whose index correspond to the rows of (x) containing the zero eigenvalues 
and, subsequently, (%) eliminating the rows and the columns containing the 
zero eigenvalues. Hereinafter (U) shall be an m x (n + 1) rectangular matrix 
with orthonormal rows and (x) an m x m diagonal and nonsingular matrix. 
We set 



w p = J2 U p,s u s, p = 1, • • • , m (35) 



n+l 

E 

s=l 
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and consider the w p as the arbitrary independent variables. Using Eq. (JT5|) 
we can write 

(Uy(x)(U) = (C) = {V)\Q){V). 

The row-spaces of (U) and (V) necessarily coincide with the (U) and (V) right 
image spaces that in turn coincide with the eigenspace of (C) associated to 
the eigenvalue zero. There exists then a non-singular mxm matrix (R) such 
that 

(V) = (R)(U). 

Now, for any complex m-tuple w = (wx,...,w m ), we have w\x)w = 
w^(RY(Q)(R)w which leads to ( X ) = (R)^(Q)(R). Thus ( X ) and (Q) are 
related by a congruence and Sylvester's inertia law 17 applies, and the num- 
ber of positive (negative) p p s coincides with the number of positive (negative) 
(x p — Xi)s. In this way the generalization of the Caratheodory theorem is 
complete. 

3 Conclusion 

Summarizing, given n complex numbers Ci, . . . , c n one considers the Hermi- 
tian Toeplitz matrix (C) defined by (j3J). One evaluates its distinct eigen- 
values, denoted by Xi < ••• < Xu with multiplicities /ii, . . . , /i„. Set- 
ting (Ci) = (C) — Xi(I) with I = l,...,v, the resulting matrices with 
/ 7^ 1 and I ^ v are indefinite. For each of these I values, the complex 
numbers — xi, c\, . . . , c n also can uniquely be written in the form (Q) with 
rrii — ( — Hi + Sg=i Hq) iff i) the rank and the principal rank of (Cj) are equal 
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[its value turns out to be equal to mi\ and ii) the (mi + 1) x (mi + 1) ma- 
trix (Si) [defined by (J2*HJ) . (fHKJ) . (jSEJ) and (|2*3|) ] is non-negative definite and 
has rank mi. In proving these results it is essential to know the conditions 
that must be obeyed by the coefficients of a polynomial equations for all its 
distinct roots to lie on the unit circle. The answer to this problem is given 
by the lemma reported in appendix A. 

As last remark we observe that, in theorems 2.2 and 3.4 of Ellis and Lay 9 , the 
assumption that the resolvent has unimodular roots can be removed by the 
aforesaid lemma. It can be substituted with the constructive requirements 
that: a) the coefficients of the resolvent equation obey Eq. ()37|) if the given 
Toeplitz matrix (T) is not Hermitian (oppositely, the condition is already 
fulfilled) , b) if the discriminant of the resolvent equation is equal to zero, one 
algebraically eliminates 18,19 all the multiple roots from the resolvent obtain- 
ing the lowest degree resolvent equation [i.e. the equation with roots equal 
to all the distinct roots of the outset resolvent], c) from the coefficients of 
the (new) resolvent equation one constructs matrix (S) defined by (|2fi|) and 
one checks its positive definiteness. In the only affirmative case the Ellis-Lay 
generalized factorization of (T) is possible. This reduces to Caratheodory's 
generalized one if (T) is Hermitian and the outset resolvent has no multiple 
roots [i.e. step b) is not required]. 
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A Unimodular roots' conditions 



We shall now prove a lemma that states the necessary and sufficient condi- 
tion for all the zeros of a polynomial equation with complex coefficients lie on 
the unit circle. The so far known theorems that ensure such property leans 
upon the existence of other polynomial with unimodular roots 10 ' 11 , while the 
following lemma only involves the coefficients of the given polynomial. 
In appendix B it is shown that, if the coefficients a\ of the iVth degree poly- 
nomial equation 

N 

P N (z)= J] Cz-ej)=E a '* l = > ( 36 ) 

l<j<N 1=0 

obey the following conditions 

aw = 1, | a.o | = 1 &nd = aAr_ m /a for m — 0,...,N, (37) 

the roots of the equation are such that ej = for j = l,...,m and 

ii, . . . , i m equal to a permutation of {1, ... , m}. The unimodularity of all the 
roots being not assured by this property, it is natural to ask: which are the 
further conditions to be obeyed by the a;s for all the roots to lie on the unit 
circle? 

The answer is given in the lemma reported later and based on Caratheodory's 
theorem. 

We first observe that it is not restrictive to assume -as we do below - that the 
roots of Eq. (J36|) are distinct because possible multiple roots can algebraically 
be eliminated 18,19 . Consider the following symmetric functions of the roots 
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of flmj) 



N 



a p = J2^ P p = 0,±l,±2,.... (38) 

5=1 

They exist for negative p integers because ao ^ 0. For non-negative ps the 
a p s are uniquely determined from the coefficients of (|3T)j) by the following 
relations (see, e.g., Ref. [18], Chap. XIII) 

a = N 

Cbj<l<J\ = — CtjV— 1 

aNO~3 +CLN-l0'2 +Q-N -2&1 = — (39) 



aAr(TAr_i+ . . . +a A a 3 +a 3 a 2 +a 2 cr 1 = —(N — l)a\ 

aN&p+N+ ■ ■ ■ +a 2 a p+2 +«lO"p+l +«0O"p =0, p = 0,1,. . . 

Owing to the condition ej = 1/e^., from (|3*8*jl follows that 

cr_ p = a p p = l,2,..., (40) 

and the last of relations (|39|) holds also true for negative p integers. In fact, 
the complex conjugate of this relation by (|37j) becomes 

Oo(o/V CTp+AT + • • • + 02 0" p+ 2 + OT Cp+l + ^0 0p) = 

fl0O"-p-7V + • • • + aAr_ 2 (T_p_2 + aAT_iCr_ p _i + aAr(X_p = 0, (41) 

and the statement is proven. The previous considerations show that all the 
<7 p s are known in terms of do, . . . , a^. 

Introduce now the (N + 1) x (N + 1) Hermitian Toeplitz matrix (S) having 
its (i, j)th element defined as 

~ rr, ;. i,j = 1,...,N+1. (42) 

Assume first that the ejs are unimodular (and distinct) and, similarly to 
(JI2J), introduce a further N x (N + 1) matrix (V) with V r , s = e/ -1 . The 
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assumed properties of the e,-s ensure that (S) = (V*)(V), that det(«S) 7^ 
and that the rank of (V) is N. The three properties in turn imply that (<S) 
is a non negative definite matrix of rank N. Then, from Caratheodory's 
theorem follows that the a p s can uniquely be written as 

N 

^ = E"W, p = 0,±l,...,±7V (43) 
3=1 

with the utjS unimodular, distinct and roots of the resolvent equation gener- 
ated by matrix (S), i.e. 



Q N (z) = Ajv _1 det 



(44) 



V 1 z ■■■ z N - x z N j 
(here An denotes the determinant of the N x N upper left principal mi- 
nor of («S)). The comparison of [)4ISj1 with (jHEJ) and the uniqueness of 
the Caratheodory decomposition imply that T\ — . . . = Tjv = 1 and 
{ui, . . . , tj n } = {e 1; . . . , €n}. From the last follows that Qjq{z) = Pn( z )- 
At this point we can state the lemma: 

the roots of an N degree polynomial equation Pn(z) = are unimodular and 
distinct iff its coefficients ai, besides obeying conditions \3T\ j, are such that 
matrix (S), defined by \4 6 A) ; ^ s non-negative definite and has rank N. 
# The necessity of the lemma has already been proven. To prove its suffi- 
ciency one has to show that the properties that the rank of (S) is iV and that 
(S) is non-negative definite ensure that the roots of Pn{z) = are distinct 
and unimodular, respectively. In fact, the first property implies that An 7^ 0. 
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From definition (jSSj) and property (j4*Uj) follows that 



N ^2 



N-U 



A N = det 



a 



(J, 



N 

Eclet 

3=1 



/ 1 



(45) 



a 



-N+3 



AT-1\ 



0", 



E det 

i<h,-,3N<N 



1 



v J 

r 



-N+l 



-N+2 



32 
N+'A 



c .il 



N-l\ 
N-2 



.IN 



'-3N 



L JJV 



1 / 



The last expression can also be written as 

f 



E 



e° e 1 e ■ 

l<3l,~;3N<N c Ji c J2 • • • c 3n 



N-l 



det 



( l e h 4 

1 e e 2 



U J1 



JV-1 • 
N-l 



\12 



V 1 9* 4 



N-l 



(46) 



Within the sum the only terms with ji ^ H ^ ■ ■ ■ ^ 3n can differ from zero. 
In other words, the possible values of {ji, . . . , jrsr} correspond to the possible 
permutations of {1, . . . , N}. The values of the corresponding determinants 
are ( — ) J=> rii<i<j<A r ( e i — e «) where P is the number of the transpositions re- 
quired for passing from . . . , j N } to {1, ... , iV}. One concludes that 



^n= n fo-*)(iA,--iAo- 

l<i<j<N 



(47) 



Thus, An 7^ ensures that the roots of Pn(z) = are distinct. We show 
now that the resolvent of (S), i.e. Eq. (JHJ), coincides with Pn(z). In fact, 



Qn(z) can be written as Qn{%) = X)p=o9pZ p = with 



■1) 



N+p 



A 



-det 



N 





• °p-l 


a P+ l 






0-1 ■ 






®N-2 ^N-l 


• (48) 






<J -N+p+2 


■ <7 a i 1 
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Manipulations similar to those performed in Eq.s (j45H46J) convert the deter- 
minant present in (|4"K|) into 

(\ ... e .P-i e .p+i ... e . N-i e .N\ 



£ ,0,1 ~ 7V^l det 



1 ... f . P-l f . P+l ... f . JV-1 , AT 

- 1 - C J2 C J2 C J2 C J2 



Vl ••• e p_1 e p+1 ••• e 7V ~ 1 e- ^ / 



Using the property that ej N = — Y,p=o a>p£j p , the above expression becomes 



— a. 



1 ... f . P-l , . P+l ... iV-l p 

^P rl„+ - 1 \72 C J2 fc j2 C J2 

.0,1 , 7V-l QeT 

i<h,-,3N<N n 32 ■ ■ ■ ; • ; • • : 

\1 ... e . p-i e . P+l ... f . Jv-i f . P/ 

and from (|46J1 and (}4*8|) one concludes that q p — a p , p — 0, . . . , N. In this way, 

the resolvent of (5) coincides with P N (z). Consequently, the e,-s also are uni- 

modular because the assumed non-negativeness of (S) and Caratheodory's 

theorem ensure that the unimodularity is true for the roots of resolvent 

Qn(z). Thus, the lemma's sufficiency is proven. # 

From the lemma follows, for instance, that the quadratic and cubic equations 
have distinct unimodular roots iff their coefficients are as follows 

N = 2: ai = pe^ /2 , a = with < p < 2 and £ [0, 2tt) 
N = 3: a 2 = pe i ^\ a 1 = pe i ^ ) a = with 
either < p < 1, 0,-0 £ [0, 2tt) 

orl<p<3, 0£[O,2tt), (20-$(p)) < 3^ < (20 + $(p)^ 
where $(p) = arccos[(p 4 + 18p 2 - 27) /8p 3 ]. 

An example of Hermitian Toeplitz matrix whose resolvent does not obey the 
conditions required by the lemma because the rank of (S) is smaller than N 
is given at the end of appendix D. 
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B Properties of the resolvent coefficients 

Given a polynomial equation of degree m 

m m 

P m ( z ) = U( z -e J ) = J2a l z l =0 

j=l 1=0 

with ao 7^ 0, we have the interesting property: 
the coefficients ai obey to 

ai = a m _i/a , l = 0,...,m (49) 

iff the roots of the equation are such that ej = for j = 1, . . . , m and 
ii, . . . ,i m equal to a permutation of {1, ... , m}. 

# This condition amounts to say that the polynomial is s elf- reciprocal 10 . To 
prove the necessity one starts from expression (|19|) of a/. Taking its complex 
conjugate and using the assumed property of the roots one finds 



m = H m ~ l E 

l<jl<-<jm-l<N 



^jl ^jm — l 



\m—l 



\m—l 



E 



e l ■ " • e m. l<j 1 <... 
( _\m— I 

E 



E 



<N e *ii ' ' ' 




ei ■ 






■ £i , 

b m — L 






e *i ' ' ' e k 





j=l tj i<j 1 <...<jj< m "0 

To prove the sufficiency one observes that 



m 

n 


i z ~ e l 


m 

) = E«7^ 


m 

= E 


i=i 










m 

t=o * 


„m m -i 

= -ne 




a 





a m-j j 



(50) 
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The previous manipulations require that no root is equal to zero and this is 
ensured by the condition ao ^ 0. With z = ej, whatever j in {1, . . . , m}, the 
first product in (150(1 vanishes. For the second to vanish one must have that 
l/e] = 6j. and the property of the roots is recovered. In passing it is noted 
that the property is true also when some roots have multiplicity greater than 
one.# 

We show now that conditions (|49jl are obeyed by the a/s defined by Eq.s (123)1 
and (J21J). In fact, setting / = in Eq.(J23J) one finds 

a = (-) m det f C 2 - " > m+1 ) /D = (-) m e ie i (51) 

where the last equality follows putting n = m + l,p=l and q = 2 in 1)59)1 . 
For the remaining / values we substitute (j2Hj) in (|4*9"J) obtaining 

det (C 1 '- ■ i i ' +2 - " - m+1 ) det (C 2 ^ - m+1 ) (52) 
= detfC ,1 '•■'• m - ^ • m - ^+2 ■"- m+1 )det(C ,1 • / = 0, • • • ,m 

where we let / take value m because ao 7^ 0. Taking n = (m + 1), 
(ji, ■■•,3m) = (1, ' • • ,m - l,m - I + 2, ■ ■ ■ ,m + 1) and (ii, ■■■ ,i m ) = 
(1, • • • , m) in (|57|) one finds that 

det f C 1 '- ."*-«.»»-'+».••■ ' m + lN ) = det fC 1 -- ' M+2 '" ' m+lN ) . (53) 

V l,-,m / V 2,— , m+1 / 

Using the property that all the mxm strictly principal minors of (C) coincide, 
the rhs of ()52|) becomes 

det f C 1 '- • u+2 "" • m+lN ) det f C 2 '"' ' m+1 ) . (54) 

V 2,-, m+1 / V 2,-. m+1/ 
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From (|5T|) follows that 



det(C 2 >~ - m+1 ) =det(\ 1 >"- m )det(C 2 - •■ m+1 

V 2,— ,m+l' V 2,— ,m+l/ V 1,— ,m 

and 

(IctfC ' " ' ' ' ) (let (A 1 '-- m )(lcl(V' : - " ^ ' ' V 

V 2,-,m+l ' V 2,-,m+l/ V l,-,m ' 

The substituion of the above two relations in (jHljl yields 



detfA l --- m 

V 2,-,m+l 



det(C7 1 '-' , ' ,+a '-' m+1 )det(C a '- ,m+1 



that coincides with the left hand side of by 

C Some properties of Hermitian Toeplitz 
matrices 

We list here a series of properties obeyed by a square Hermitian Toeplitz 
matrix (C) of order n and partly reported in Ref. [8]. 

(a) - Its elements obey to 

C r ^ s Cs,r Cr—s Cs—ri S 1, . . . , (*-^) 

so that all the elements of (C) contained in a line parallel to the main diagonal 
are equal. 

(b) - One has the reflection symmetry with respect to the second diagonal 
formalized by the condition 

CV,s = c n +i_ Sjn -|_i_ r . (56) 

(c) - All the (m x m) principal minors of (C), whatever the considered rows 
(and columns), are identical. 
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# The property is a consequence of (a).# 

(d) - For any choice of m rows (1 < i\ < ■ • ■ < i m < n) and m columns 
(1 < ji < ■ ■ ■ < j m < n) it results 

detfC^-'M = detfC T <™+ 1 -^),-,(™+i-n) ) = detfc ("+!-*»),- ,c»+wi) \ 

V H,-,W V (n+l-i m ),- ,(n+l-hy \ (n+l-i m ),- ,(n+l-ii)^ 

(57) 

# The first equality, where (C T ) denotes the transposed of (C), follows from 
property (b) and the second from the Hermiticity of (C). # 

The following properties, that we think to be original, hold only true for 
Hermitian Toeplitz matrices having their rank equal to the principal one. 

(e) - If the principal rank of a Hermitian Toeplitz matrix (C) is equal to 
the rank m(< n) of (C), any (m x m) strictly principal minor of (C) is non- 
singular. 

# The property immediately follows from the definition of "principal rank" 
of a matrix, reported above Eq. (JHJ), and (c)^. 

(f ) - For any (n x n) Hermitian Toeplitz matrix of rank equal to its principal 
rank m(< n), the determinant of any of its minors formed by m subsequent 
rows and m subsequent columns is simply related by a phase factor to the 
determinant of the (strictly) principal minor contained in the considered rows 
or columns, i.e. 

det (C {q+1) '- ' q+m ) = e i0 ^ q det (c p+1 '- ' p+m ) , (58) 
det (C (9+1) '- ' q+m ) = e idp -* det (C ' 9+m ) , 

V p+l,--- ,p+m/ \ ty+l,--- ,g+m/ 

# P - 9 e K, p,q = 0, 1, . . . ,n - m. 
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# Clearly, if the first of the above two equalities is true the second also 
is true because of (c). To prove the first of equalities (|HH)l one observes 
that (d) implies that any m distinct rows of (C) can be written as linear 
combinations of m other distinct rows (see, e.g., Ref. [18], Chap. III). Hence 
rows (p+1), . . . , (p+m) can be expressend in terms of rows (q+1), . . . , (q+m) 



as 



q+m 

C r ,s= K,tC t ,s, r = (p + 1), . . . , (p + m), s = l,...,n, (59) 

t=q+l 

where the X Ttt are suitable numerical coefficients. From these relations follows 
that 

det (C p+1 '- • p+m ) = det f A ' 9+m ) det (c p+1 <- > p+m ) . (60) 

Due to (e) the left hand side of fl6()(l is different from zero so that both factors 
on the rhs are different from zero. The complex conjugation of fl60j) . by the 
Hermiticity of (C), yields 



det(C p+1 '-' p+m = det(A 9+1 '-'« +m det C 9+1 '"' 9+m . (61) 

V p+l,---,p+m/ V p+l,---,p+m/ V p+1,--- ,p+m' 



From Eq. (|59J) also follows that 



detfC 9+1 '''' 9 + m ) = detfA 9+1 '-'« +ro )detfC 9+1 '-' 

V p+l,---,p+m/ V p+l,---,p+m/ V q+1,-.. , 



q+m 

q+m 



= det (A ) det(C p+1 ' 

V p+1,— , p+nv V p+1,— , p+m 

where the last equality follows from (c). The substituion of the last equality 
in Eq. (|50|l and the fact that, by assumption, det(C p+1 '"'' p+m J ^ imply 

» p+1,--- ,p+m/ 

that 

s 2 

= 1 p, q — 0, . . . , n — m, (62) 



det (A q+h -' q+m 

V p+1,-- , p+m 
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and Eq. ()58|) is proven. That the phase factor depends on p — q instead of 
(p, q) follows from the fact that the two determinants present in do not 
change with the two substitutions p — > p + 1 and q — > q + 1 owing to (JH%|) .# 
An immediate consequence of (f) is the property that 

(g) - any (m x m) minor formed by m subsequent rows and m subsequent 
columns of a Hermitian Toeplitz matrix with rank equal to its principal rank 
m is non-singular. 

D Numerical examples 

To illustrate the application of the results discussed above we report three 
numerical examples. 

1 - The first shows where it is impossible to express a set of CjS in 

terms of positive and negative pjS. Assume that c\ — 0, C2 = 0, C3 = 1. 
The corresponding matrix (C) has eigenvalues equal to —1, 0, 0, 1. The 
matrix (Ci), obtained by setting Co = 1, is semi-positive definite with rank 3 
and eigenvalues equal to 0, 1, 1, 2. The solutions are: pi = p 2 = p% = 1/3 
and ei = 1, 62 = — e i7r//3 , 63 = e* 27r//3 . The matrix (C2), obtained by setting 
Co = —1, is semi-negative definite with rank equal to 3 and eigenvalues 
equal to —2, —1, —1, 0. The solutions are: p\ = p 2 = P3 = _ 1/3 and 
ei = —1, e 2 = e" 1 "/ 3 , e 3 = — e l2n ^ 3 . Finally, the matrix (C 3 ) obtained by 
setting Co = coincides with (C). It is non-definite, has rank equal to 2 and 
principal rank equal to 0. For this reason it is impossible to write 0, 0, and 

1 in the form (1) with m = 2 as it is easily checked. 

2 - The second example considers the case where c\ — 1, c 2 = 0, c 3 = 1. The 

26 



eigenvalues of the associated matrix (C) are —2, 0, 0, 2. Setting c = the 
resulting (C) matrix coincides with (C). It is non-definite, its rank is two and 
equal to its principal rank value. The generalized Caratheodory's theorem 
applies and the solution is e\ = 1, e 2 = — 1, pi = 1/2 and p 2 = —1/2. With 
these values one easily checks that c = p\ + P2 = 0, C\ = pi€i + p 2 e 2 = 1j 
c 2 = Pie? + p24 = and c 3 = Pie? + P2ei = 1. 

Setting Co = 2 the resulting (C) is non-negative defined with rank 3. The 
solution is: e\ = 1, e 2 = i, €3 = —i, p\ = 1, p 2 = 1/2 and P3 = 1/2. The last 
choice Co = — 2 defines a non-positive defined (C) matrix with rank equal to 3 
with solution: ei = —1, e 2 = z, €3 = — i, pi = —1, p 2 = —1/2 and P3 = —1/2. 
3 - The last example corresponds to have Co = So, c\ = 5o + iS\t %Lp and 
c 2 = (So + 2z5i)e 2ll/ ' with So, Si and <p reals. The eigenvalues are and 
(3#o/2)(l ± ^/l + 85f /(3<5q))- The Hermitian matrix (C) is indefinite and has 
rank 2 under the assumption that no further eigenvalue is equal to zero. The 
resolvent equation is P(z) = z 2 — 2e lip z + e 2lip = so that its coefficients 
ao = e 2lip , a\ = —2e ttp and a 2 = 1 obey Eq. (|3*7|) . From (39) one immediately 
finds that o"o — 2, o\ = 2e lip and cr 2 = 2e 2lip and the resulting matrix (S) has 
rank 1. Hence the previous resolvent equation has a unimodular root with 
multiplicity 2: in fact, P(z) = (z — e tlf ) 2 = 0. According to our analysis, 
(C) cannot be written as (V)(A)(V)^ with (A) equal to a diagonal matrix. 
In this example, however, the roots are unimodular and (C) factorizes in the 
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Ellis-Lay form 9 as 
/ 



So, 

2iS 1 )e- 2itfi , 



(8 + iSJe*, 
So, 

(8 - iSJe-*, 



(S + 2iS 1 )e 2i ^ 
(So + iSJe* 
So 
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